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The long-range non-additive three-body dispersion interaction coefficients Zm, Z\\2, ^113, 
and Z122 are computed for many atomic combinations using standard expressions. The atoms 
considered include hydrogen, the rare gases, the alkali atoms (up to Rb) and the alkaline- 
earth atoms (up to Sr). The term Zm, arising from three mutual dipole interactions is 
known as the Axilrod-Teller-Muto coefficient or the DDD (dipole-dipole-dipole) coefficient. 
Similarly, the terms Z112, Z113, and Z122 arise from the mutual combinations of dipole (1), 
quadrupole (2), and octupole (3) interactions between atoms and they are sometimes known, 
respectively, as DDQ, DDO, and DQQ coefficients. Results for the four Z coefficients are 
given for the homonuclear trimers, for the trimers involving two like-rare-gas atoms, and for 
the trimers with all combinations of the H, He, Li atoms. An exhaustive compilation of all 
coefficients between all possible atomic combinations is presented as supplementary data. 

PACS numbers: 31.15.ac, 31.15.ap, 34.20.Cf 



The DDD (dipole-dipole-dipole) or Axilrod-Teller-Muto interaction [1( is a non-additive dis- 
persion interaction that arises between three polarizable atoms or molecules. In the long-range 
potential energy function it enters at the order of the inverse ninth power of the characteristic 
separation of the three particles. It is studied in the physics and chemistry of gases, liquids, 
adsorption of atoms at surfaces, and materials science, and more recently for applications to ultra- 
cold atoms. In addition to the DDD interaction, there are terms arising from mutual interactions 



I. 



INTRODUCTION 



of higher-than-dipole induced multipoles , such as the dipole-dipole-quadrupole (DDQ), the 
dipole-dipole-octupole (DDO), and the dipole-quadrupole-quadrupole (DQQ) interactions. Values 
of the coefficients are useful for creating models of potential energy surfaces and for comparison 
with more sophisticated methods that include exchange energies. 

For a variety of reasons, there has been increased interest in the non-additive three-body dis- 
persion coefficients. Some of the primary motivations come from studies of ultra-cold atoms. For 
example photo- and magneto-association techniques are being developed to form homonuclear and 
heteronuclear trimers. The collisions that govern such events typically occur at low energies where 
the importance of dispersion events are enhanced and there have been a number of studies that 
have investigated the impact of three-body non-additive interactions ^~2\- Another application 
area is in fine-tuning the interaction potentials of rare gase s for the description of their bulk ther- 
mophysical properties such as the virial coefficients |s-12] and for computing properties of their 



solid forms 
in rare gases 



13K There has also been interest in describing the structures of metal atoms embedded 
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18] . Finally, the dispersion interactions give reference data that could be utilized 



to validate ab initio calculations of the full three-body potential surfaces at large inter-nuclear 
separations [191421 1. 



The DDD coefficients are known very accurately for systems involving combinations of H, He and 
Li through ab initio calculations 



22 



261 ] . while the most accurate DDD values for the heavier rare 



gases are derived from pseudo-oscillator strength distributions (POSDs) u sing 



and other experimental data 



27 



291 ] . Doran and Zucker [3j] and Doran 3CJ, [3l| calculated higher 



dipole polarizabilities 



order coefficients for the rare gases using Hartree-Fock methods. Their results were surpassed in 
accuracy by those of Thakkar et al |32l] . Reasonably accurate values for the DDD coefficients were 

3411 and for the alkaline-earth metal atoms 34J]. An alternate 



computed for the alkali-metal 
)pr 

M 



app roach to the determination of three-body dispersion coefficients is to use combination rules 
35H38I]. which utilize estimates of other properties such as polarizabilities and two-body dispersion 



coefficients. 

The present work represents the most comprehensive evaluation of three-body dispersion co- 
efficients so far presented in the literature. Coefficients are computed using oscillator strength 
sum rules up to a combined inverse power of 15 in the interspecies separation distance are given. 
Extremely precise dispersion coefficients are given for all combinations of the H, He and Li ground 
states since definitive results for these atoms can often serve as benchmarks to test methodologies 
used for more complicated systems. In addition, dispersion coefficients are given for all the three- 
body homonuclear trimers involving the rare gases, the alkali atoms up to Rb, and the alkaline 



3 



earth atoms up to Sr. Space requirements preclude the presentation of data for all possible trimers 
involving heteronuclear systems. Tabulated coefficients are only given for systems with two like 
rare gas atoms. The three-body dispersion coefficients for almost all other possible three-body 
systems are detailed in the supplementary data submitted with this paper. 



II. DEFINITIONS OF THE DISPERSION COEFFICIENTS 
A. Expressions using oscillator strength sum rules 

The leading non-additive dipole-dipole-dipole dispersion interaction between three atoms (a, b, 
c) can be written as 

Viii (R a , Rb , R c ) = Zfii Wiii (R a , Rb , R c ) , ( 1 ) 

where Zfii is the three-body dispersion coefficient and the function Wf^Ra, Rb> Rc) is detailed 
in [l, j]], e.g. 

u/aW-R n r ^ 1 + 3cos(fl ab ) cos(gbc) cos(fl ca ) 

Win (Ka,Kb,KcJ - 3 R 3 r3 ' W 

where # a b is the interior angle at c due to atoms a and b, i? a b is the distance between atoms a and 
b, and 

f (i) f (i) f (i) 



yabc _ ° 
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The sums over i,j, and k involve the dipole oscillator strengths f^' and the energy difference eoi 
from the ground state of each atom. There are alternate definitions of Zm where Zm is a factor 



of 3 smaller with the Wm geometric factors having an extra multiplying factor of 3 % 
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32]. 



More generally, the three-body potential for arbitrary multipole combinations can be written 

V tjklc ( Ra ' Rb ' Rc ) = Z iaVc W ti b ic ( R a> R b, Rc) , (4) 



where the general dispersion coefficient is written 

So) Mb) Ate) 
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(5) 
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The /o* refer to t 



nc mu 



the definitions in 2f], [34], e.g. 



tipole oscillator strength from the ground state of each atom according to 



(21 + 



(6) 



The next term in the three-body dispersion series is the dipole-dipole-quadrupole (DDQ) con- 
tribution. In the case of a heteronuclear trimer there are three separate terms since the quadrupole 
excitation can be associated with any of the atoms, a, b, or c. One term is 



^ml^a, Rb, R-c) — ^mlR-a, Rb, Rc) • 



abcyrrabc i 



(7) 



The Zf& coefficient is 



L \\2 



f (l) f (l) f (2) 
■/a,0i-/b,0jJc,0fc 

2 ~ ea,0ieb,0jr'ec,0A; 



IE 



(£a,0i + £b,0j + e c,0fc) 



(£a,0i + £b,0j)(ea,0i + €c,0fc) (eb,0j + ^cflk) 

The spatial factor (dropping the arguments for compactness) for Wf^ 1S 



(8) 



Ti/abc 
^112 



9cos(0 ab ) - 25cos(3# ab ) 



6 COs(6> bc -#ca) [3 + 5cOs(26> ab )] 



0) 

There is no need to give the Z coefficients for the other combinations. The spatial factors for VFf^ 

as 
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Wfi§ and iy 2 a 2 b 2 c can be written j^fl] 



Tirabc 
^122 
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3cos(6*bc) + 15cos(3# bc ) 



rirabc 
"113 



+ 20cos(6> ca -# ab )[l -3cos(20 bc )] 
+ 70cos[2(6> ca -# ab )] cos(# bc 
5 
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+ 6 cos(6 l bc - 6> ca ) [9 cos 6» ab + 7 cos(36» ab ) 



(10) 



(11) 



Wi 



abc 



222 



USR^RtR^ 



490 cos(26> ab ) cos(26> bc ) cos(26» ca ) 



- 27 + 220 cos(0 ab ) cos(# bc ) cos(0 ca ) 

+ 175{ cos[2(0 bc -fl ca )] + cos[2(0 ca -0 ab )] 



+ cos[2(0 ab -0 bc )] 



(12) 
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TABLE I: The dispersion constants (in a.u.) for homonuclear trimers consisting of either hydrogen, helium 
or lithium. All other dispersion coefficients are given by the symmetry conditions as described in the text. 
The numbers in parentheses are the computational uncertainties arising from incomplete convergence of the 
basis set. Values for hydrogen are correct to all quoted digits. The numbers in the square brackets denote 
powers of ten. 

System Zm Zm Z122 Zm ^222 

H 21.642464510636 78.707166385334 289.37119056926 712.03297892387 1077.4737895329 

He 1.47955860643(1) 2.77280376179(1) 5.21831397022(6) 12.37878701399(5) 9.8654795914(2) 

Li 1.70615(2)[5] 1.78338(2)[6] 1.99321(1)[7] 5.12903(4)[7] 2.463549(6) [8] 



Expressions for rearrangements of the different multipole excitations amongst the different atoms 
can be written down transparently. The expressions for the homonuclear trimers, i.e. and 
Wffi* do not require any modifications to the formulae listed above. 

B. Combination rules 

The present approach to the calculation of the dispersion parameters is based on explicit eval- 
uation of oscillator strength sum rules. An alternative approach is to use integrals over the po- 
larizabilities at imaginary frequencies [jj, [3] . Approximate expressions that use existing calcula- 
tions of polarizabilities and two-body dispersion parameters in combination rules have been derived 



35, 



361 ] . One of these combination rules is known as the Midzuno-Kihara approximation 35j, |36J] . 



For a trimer consisting of three identical atoms, one has 

C 9 

Other more complicated ex pre ssions exist for heteronuclear trimers 



(13) 



rules have been developed 



35 



361 ] and other combination 



381 ]. Estimates using the Midzuno-Kihara approximation for Zm are 



given in some tables and abbreviated as MK. 



III. UNDERLYING DESCRIPTION OF ATOMIC STRUCTURES 



The evaluation of the three-body dispersion parameters involves the construction of tables of 
oscillator strengths to describe the excitation spectrum of all atoms considered. These oscillator 
strength distributions were constructed by a variety of means depending on the specifics of each 



6 



atom. A short description is now given for each atom since the accuracy of the dispersion coefficients 
depends on the details of the underlying structure model. 

A. Hydrogen 

The hydrogen atom and its excitation spectrum was computed by diagonalising the hydrogen 



hamiltonian in a basis of Laguerre type orbitals (LTOs) 39|,|40]. Such a basis has the property that 
it can be expanded to completeness without any linear dependence issues. Very accurate dispersion 
parameters can be obtained with relatively small basis sets. 

B. Hylleraas descriptions for He and Li 

The helium atom wave function is expanded as a linear combination of Hylleraas functions. 
The resulting oscillator strength sum rules are capable of giving the dipole polar izability to an 
accuracy of 12 significant figures. The lithium atom wave function can also be represented as a 
linear combination of Hylleraas functions. The achievable accuracy for the Li ground state dipole 
polarizability is 5 significant figures. Details of the Hylleraas method and its application to He and 
Li can be found in Refs 



13, 



26]. 



C. Pseudo-oscillator strength distributions for heavier rare gases 

The pseudo-oscillator strength distributions (POSDs) for the heavier rare gases, Ne, Ar, Kr 
and Xe, come from two sources. The dipole distributions are taken from Kumar and Meath 



27|, 



2g |, which were constructed by reference to polarizability data and other experimental data. 
The quadrupole and octupole pseudo-oscillator strengths are based on Hartree-Fock (HF) (or 
relativistic HF) single particle energies and radial expectation values. These were given tunings to 
reproduce calculations of higher multipole polarizabilities and homonuclear dispersion coefficients 



4l| 



D. Semi-empirical description of structure for alkali and alkaline-earth atoms 



The transition arrays for the alkali atoms are those which were used in calculations of the 
dispersion interactions between these atoms and the ground states of hydrogen and helium 42|. 
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TABLE II: The three atom dispersion constants for combinations of H, He and Li with two like atoms. The 
other dispersion constants are given by the symmetry relations Z^fi—Zf^i , ^212— ^122 > an d ■ZJfi = ^ifi ( m 
a.u.). The numbers in parentheses are the computational uncertainties arising from the finite basis set size. 



System 


■^111 


ZI12 




^122 


H-H-Hc 


8.10224087465(2) 


14.6777067844(1) 


30.3021391145(1) 


55.1425238705(2) 


H-H-Li 


275.987(3) 


4195.023(6) 


944.97(1) 


14670.29(2) 


He-He-H 


3.26806489632(2) 


12.7469884056(1) 


5.99712674013(3) 


23.5605644280(1) 


He-He-Li 


29.8259(6) 


504.658(1) 


53.228(1) 


907.467(2) 


Li-Li-H 


6133.7(1) 


20621.1(2) 


P^£"* /""v A C~\ C\ { /~*\ 

72048.9(6) 


243884(2) 


Li-Li-Hc 


1917.48(2) 


3394.56(4) 


22845.5(2) 


40500.9(4) 


System 


yaab 
^221 


C^aab 
^113 


^aab 
^131 


7aab 
^222 


H-H-Hc 


115.1382084748(3) 


64.4548538033(2) 


277.5076627177(7) 


210.661799531(1) 


H-H-Li 


3245.34(5) 


135490.15(9) 


8349.17259(1) 


51577.89(8) 


He-He-H 


11.0328224804(1) 


118.9339509877(6) 


26.4993528199(1) 


43.6808551603(3) 


He-He-Li 


95.057(2) 


17139.665(9) 


232.194(6) 


1633.531(4) 


Li-Li-H 


966692(2) 


180949(2) 


214091(2) 


3312094(6) 


Li-Li-Hc 


314438.8(7) 


14753.5(2) 


682925(7) 


558901(2) 



The transition arrays for the alkaline-earth atoms are those which were used in calculations of the 



3|. 



dispersion interactions involving the alkali- and alkaline-earth atoms 

The wave function is partitioned into core and valence parts, where the core is described by a 
HF wave function. The ground and excited state pseudo-spectrum are computed by diagonalizing 
the fixed-core Hamiltonian of one and two electron functions written as a linear combination of 
LTOs. The HF core Hamiltonian is supplemented with a semi-empirical core polarization potential 
tuned to reproduce the energies of the low-lying spectrum. For all practical purposes the results 
are converged with respect to increasing the dimension of the orbital basis. 

Excitations from the core are included in the calculations of the dispersion coefficients by the 
construction of pseudo-oscillator strength distributions based on HF energies and expectation val- 
ues 3j] . This approach to the determination of atomic structure is referred to as the configuration 



interaction plus core polarization (CICP) method in the remainder of this article. 
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IV. NUMERICAL RESULTS 
A. H, He and Li 

A complete set of dispersion coefficients for all combinations of the light systems involving H, He 
and Li up to a total inverse power of R 15 are given in Tables HI ITT1 and ITTT1 The dispersion coefficients 
for these systems are given special prominence since they can be used to validate high precision 
quantum chemistry calculations of potential energy surfaces for these trimers. Table IIVI compares 
values of Zm with CICP results to give an accuracy benchmark for the CICP calculations. 

Certain symmetry conditions apply when the trimer consists of three homonuclear atoms. One 
can write Zfff =Zfff =Z?ff , Zfff =Zfff =Zgf , and Zfff =Zgff =Zff§. The results for the H, He and 
Li homonuclear trimers are given in Table HI The dispersion coefficients for hydrogen are given to 
14 digits. This does not represent the achievable precision. It is possible to compute the Z^ 2 i a to 
30 significant digits by using quadruple precision arithmetic and a basis of 30 Laguerre functions. 



The Zi ± i 2 i 3 coefficients are in perfect agreement with those of Koga 25|. The Koga values use 
different conventions for relating the Z^ 2 i 3 to the geometric functions and are not quoted in Table 

El 

For helium, a precision of at least 10 significant digits was achieved, while for lithium six digits 
were tabulated with some uncertainty usually occurring in the sixth digit. There have been no 
previous high precision calculations of these coefficients except for Zm which were previously 
computed by Yan et al using the Hylleraas basis {^(J . The present coefficients are in harmony with 
these earlier calculations. 

The uncertainties in the dispersion parameters were estimated by an examination of the conver- 
gence for a series of successively larger calculations 



17 



26] . Relativistic and finite mass corrections 



for these atoms can be expected to start to occur by the fourth to sixth digits. Therefore the pre- 
cision achieved by the present non-relativistic calculations is sufficient for any practical calculation 
since finite-mass and relativistic corrections would need to be incorporated to further improve the 
overall level of accuracy. 

The Zm coefficients for all H, He, Li combinations are given in Table HVl They are compared 
with the earlier Hylleraas calculations of Yan et al 261 and coefficients from from CICP type cal- 

n 

culations [341 ] . The two sets of Hylleraas for trimers containing just H and He are generally in 
agreement when the quoted uncertainties are taken into consideration. However, the present dis- 
persion coefficients for helium are about one order of magnitude more precise than the calculations 
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TABLE III: The dispersion coefficients for the heteronuclear H- He-Li trimer (in a.u.). The numbers in 
parentheses are the computational uncertainties arising from incomplete convergence of the basis set. 





Values 




Values 




Values 


Z\w 


89.833(1) 


Z\\2 


1430.195(3) 


Z\2l 


159.785(3) 


Z211 


309.922(6) 


Z\22 


2557.314(4) 


Z212 


5068.475(9) 


Z221 


551.89(1) 




47174.62(3) 


Ziai 


695.90(1) 


Z311 


2747.29(5) 


Z222 


9081.12(2) 







TABLE IV: Comparison of the Z\\\ (in a.u.) parameter for all combinations of the trimers formed by the H, 
He, and Li with some earlier calculations. The numbers in parentheses are the computational uncertainties 
arising from incomplete convergence of the basis set. 





Hylleraas 




CICP 


System 


Present 


Ref. [26] 


Ref. [34] 


H-H-H 


21.64246451063597 


21.642464510635 


21.6425 


H-H-Hc 


8.10224087465(2) 


8.1022408743(2) 


8.103 


H-H-Li 


275.987(3) 


275.979(7) 


276.0 


He-He-He 


1.47955860643(1) 


1.4795586063(1) 


1.4798 


He-He-H 


3.26806489632(2) 


3.2680648961(1) 


3.269 


He-He-Li 


29.8259(6) 


29.824(5) 


29.83 


Li-Li-Li 


170615(2) 


170595(6) 


170873 


Li-Li-H 


6133.7(1) 


6133.5(5) 


6139 


Li-Li-Hc 


1917.48(2) 


1917.27(5) 


1919 


Li-He-H 


89.833(1) 


89.830(5) 


89.85 



by Yan et ctl 



26] 



In a few instances involving the Li atom the present Z\\\ values lie outside the error bars of the 
earlier Hylleraas calculation. The basis sets used in the Yan et al calculation were much smaller 
than those used herein. Agreement between the CICP and Hylleraas calculations is good, with no 
differences exceeding 0.2%. The older Z\\\ values of 0, l^j] and 22] are not listed since they are 
not expected to have the precision of the present calculations. 

The relative importance of the Zm coefficient with respect to the two-body dispersion coeffi- 
cients is easy to estimate and the Li-Li-Li trimer is used as an example. The two body C» and C%q 
dispersion coefficients are C 8 = 8.3429(1) x 10 4 a.u. and C w = 7.3725(2) x 10 6 a.u. [lj]. When 
three Li atoms form an equilateral triangle the spatial factor of Eq. ([2]) is ^-i? -9 . At R = 30 ao 
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the ratios of Z\x\W\\\ to C$R~ 8 and Cio-R -10 are 0.28 and 0.95 respectively. The nonadditive 
three-body interaction should be included in high-precision analysis of photo-association spectra 
for homonuclear alkali-metal trimers. These ratios are 0.0096 and 0.1661 respectively for a helium 
trimer in an equilateral triangle configuration at R = 15 ao- 



TABLE V: The three-body dispersion coefficients, Zm, ^112, -^122 and .Z113 (in atomic units) for homonu- 
clear trimers. The /-value distributions for H, He and Li use Laguerre type orbitals (H) or Hylleraas basis 
functions (He and Li) to describe the ground and excited state spectra. Dispersion coefficients for the H, 
He and Li trimers are given to additional significant digits in Table HI The heavier gas /-value distributions 
use pseudo-oscillator strength distributions 



27 



28, 



44j . while those for the other atoms come from CICP 



35 



361 ] use the present oscillator strength 



calculations. Results in the Midzuno-Kihara (MK) approximation 
distributions to compute the underlying ay and C@ needed as input. The numbers in the square brackets 
denote powers of ten. 



Systems 




Zm 








-Z112 


Z\22 


^113 


Present 


Other 




MK 


Present 


Other 


He 


1.4796 


1.47956 


a 


1.516 


2.77280 


2.77280 b 


5.2218 


12.378 


Ne 


11.95 


12.37 c 




12.78 


33.73 


34.20 d 


95.21 


189.7 


Ar 


518.3 


521.7 c 




534.4 


2.543[3] 


2.554[3] a 


1.251 [4] 


2.535 [4] 


Kr 


1.572 [3] 


1.697[3] 


c 


1.631 [3] 


9.556[3] 


10.25 [3] c 


5.827[4] 


1.174[5] 


Xe 


5.573[3] 


6.246 [3] 


c 


5.822[3] 


4.582 [4] 


5.163[4] c 


3. 778 [5] 


7.298[5] 


H 


21.624 


21.642 


b 


21.93 


78.707 


78.72 6 


289.37 


712.03 


Li 


1.7062 [5] 


1.701 [5] 


d 


1.717[5] 


1.786 [6] 


1.750[6] d 


1.997[7] 


5.138[7] 


Na 


1.892[5] 


1.906[5] 


d 


1.908[5] 


2.544[6] 


2.353[6] 9 


3.553[7] 


7.844[7] 


K 


8.318[5] 


8.493 [5] 


d 


8.488 [5] 


1.645[7] 


1.493 [7] 9 


3.354[8] 


6.212[7] 


Rb 


1.063 [6] 


1.097[5] 


d 


1.097[6] 


2.455[8] 


2.230[8] 9 


5.678[8] 


9.704[8] 


Be 


5.973 [3] 


6.135[3] 


e 


6.023 [3] 


5.335[5] 


5.535[4] e 


4.865[5] 


7.420[5] 


Mg 


3.338 [4] 


3.36 [4] 


e 


3.369[3] 


4.182[5] 


4.26[5] e 


5.307[6] 


7.667[6] 


Ca 


2.557[5] 


3.25[5] 


e 


2.626[5] 


5.131[6] 


5.34[6] e 


1.034 [8] 


1.125[8] 


Sr 


4.753[5] 


4.903[5] 


f 


4.903 [5] 


1.111[7] 




2.601 [8] 


2.983[8] 


"Independent Hylleraas 


calculation 2C 


1 













6 Pseudo oscillator strength distribution basednon accurate oscillator strength sum rules [4: 
c Many body pertuAation theory (MBPT) [331 
^Model potential [33| 

e Pseudo pscillator strength distributions (POSD). Value given is mean of the lower and upper limits. 
'MBPT |31[ 

9 Asymptotic wave function 



11 



B. Homonuclear trimers 



Dispersion coefficients for the homonuclear trimers are given in Table IVl The atoms presented 
are hydrogen, the rare gases, the alkali atoms, and the alkaline earth atoms. Values of Zm for 



22 



23, 



26, 



29, 



31 



34, 



45|448|] 



most of these atoms have been published previously 

The set of values for the rare gases heavier than He should be regarded as the current benchmark. 
The Zm coefficients use the pseudo oscillator strength distributions of Kumar and Meath 271. |28| 
and should be accurate to about ±1%. The best ab initio calculation is the MBPT calculation 
32 1 which should be regarded as superseding earlier calculations Q, 



31 



,1491, 



1 



The MBPT 



calculation overestimates the dipole polarizabilities of Kr and Xe and this leads to the MBPT Zm 
values being too large for these atoms. 

The present Z\\i values are slightly smaller than the MBPT values. The present values should be 
adopted as the recommended values for reasons outlined above since the most important multipole 
excitation is still the dipole and the pseudo oscillator strength distribution of Kumar and Meath 
is still the preferred value. The MBPT calculation 32] also gave values for Zyii and Zm,. These 



are not tabulated for reasons of brevity since they repeat the same pattern seen for the Zm and 
Z112 coefficients. The MBPT are about 5-10% larger than the present values. 

A detailed comparison of the CICP alkali atom Zm coefficients with those of a model potential 
calculation 



331 ] listed in Table [V] has been made previous 



The CICP Zm coefficients were 



331 ] did not properly incorporate the con- 



to be preferred since the model potential calculation 
tributions from the core. Two-body dispersion coefficients from the CICP calculations reproduce 
ligh quality relativistic many-body perturbation theory estimates to a high degree of precision 



34, 



51 



54j. 



For the alkali atoms there have also been estimates of the three-body dispersion coefficients 
based on a simple model pot ential tuned to give the correct binding energies for the lowest state in 



each angular momentum 



45|, 



481]. The Z112 results of Huang and Sun [45|] are listed in Table [V] The 



Huang and Sun Zm tend to be 5-10% smaller than the present CICP values. Taking the Na-Na-Na 
trimer as a specific case, the Huang and Sun Z U2 of 2.353 x 10 6 a.u. value is 7% smaller than the 
present value. A similar situation also occurs for the two-body C% coefficient for the Na-Na dimer. 
The CICP calculations gives 1.159 x 10 5 a.u. 



34], 



551 ] while the Huang and Sun calculation gives 



1.083 x 10 5 a.u. There are similar discrepancies with the Patil and Tang [48j tabulation. Patil and 
Tang give 2.35 x 10 6 a.u. for the Na-Ne-Ne Zm, which is 7% smaller than the present value. The 
present CICP calculations expose the limitations of the structure models underlying the Huang 



12 



and Sun [45] and Patil and Tang 



48] tabulations. The Patil and Tang three-body coefficients, Zm 



and were used as input in an attempt to make a global triatomic potential surface for three 
lithium atoms [5]. 

The present CICP Zi x i^ z dispersion coefficients should be regarded as the reference values. For 



23| listed in Tables ED are 



example, the Zm values from the Standard and Certain compilation 
taken as the mean values of the lower and upper bounds. The oscillator strength distributions 
for some of the atoms in the Standard and Certain compilation rely on much older and smaller 
calculations. 



C. Impurity atoms embedded in the rare gases 

The three-body dispersion coefficients for the X-He-He systems given in Table IVH where X 
denotes H, the alkali atoms up to Rb and the alkaline earth atoms up to Sr. The quadrupole 
excitation occurs on one of the He atoms for the Z^ cHc = Z™ e ^ e coefficients. The quadrupole 
excitation is associated with the impurity X atom for the Z^ eHe coefficient. 

The three-body dispersion coefficients for the X-Ne-Ne and X-Ar-Ar systems are given in Table 
IVIII while those for the X-Kr-Kr and X-Xe-Xe systems are given in Table IVIIII The quadrupole 
excitation is associated with the impurity X atom for the Z^ gRg coefficient. 



D. Other systems 



The sheer number of possible heteronuclear three-body dispersion coefficients precludes their 
listing and discussion in the present manuscript. The three-body dispersion coefficients for every 
possible three-body trimer formed from the atoms studied in this paper (H, He, Ne, Ar, Kr, Xe, 
Li, Na, K, Rb, Be, Mg, Ca, Sr) are given in the supplementary data published with this paper. 
Previous compilations of the Z^f three-body coefficients have been given for all possible alkali- 



atom trimers from Li to Cs 



33 



34 1 . In addition, Zf^f coefficients for all possible alkaline-earth 



atoms from Be to Sr have also been published 34|. The supplementary data go well beyond these 
existing tabulations in that the higher-order Zf^j ^113 j an d ^122 are a l so gi ven - 



E. The Midzuno-Kihara approximation 



The Midzuno-Kihara (MK) approximation 35], [36] for Zm has been given for all combinations 
in Tables [V] IVIIIVIII and IVIIII The coefficients in the MK approximation use the present oscillator 



13 



TABLE VI: The three-body dispersion coefficients (in atomic units), for the X-He-He systems containing 
two helium atoms. The sources for the /-values distributions are the same as those in Table [V] The 
dispersion coefficients involving H and Li are given with additional significant digits in Table [TTJ Results 
in the Midzuno-Kihara (MK) approximation 



use the present oscillator strength distributions to 
compute the underlying ay and Cq needed as input. The numbers in the square brackets denote powers of 
ten. 



Systems 


■Z111 

Present 


MK 


^112 


•^211 




■^221 


-^113 


•Z311 


Ne 


2.955 


3.067 


5.561 


8.335 


10.51 


15.69 


24.86 


46.96 


Ar 


10.25 


10.52 


19.04 


50.99 


35.51 


91.11 


84.65 


505.1 


Kr 


14.62 


15.05 


27.09 


90.66 


50.37 


168.6 


120.3 


1099 


Xe 


21.74 


22.46 


40.15 


183.2 


74.35 


339.4 


177.9 


2884 


H 


3.269 


3.324 


5.999 


12.75 


11.04 


23.57 


26.52 


118.9 


Li 


29.83 


30.13 


53.25 


505.1 


95.11 


908.5 


232.6 


1.716[4] 


Na. 


33.52 


34.30 


59.99 


624.8 


107.5 


1123 


262.4 


2.279[4] 


K 


50.73 


52.82 


90.97 


1276 


163.4 


2287 


398.1 


6.014[4] 


Rb 


56.80 


59.75 


106.0 


1538 


183.6 


2753 


447.0 


7.676 [4] 


Be 


16.54 


16.77 


29.90 


176.3 


54.13 


321.2 


131.3 


2749 


Mg 


27.05 


27.56 


48.78 


392.2 


88.07 


710.1 


213.9 


8249 


Ca 


46.34 


47.87 


83.42 


1008 


150.4 


1814 


365.7 


3.047[4] 


Sr 


55.97 


58.23 


100.8 


1354 


181.8 


2433 


442.0 


4.629[4] 



strength distributions to compute the underlying and C% needed as input. In every instance 
the MK results are larger by an amount between 1-4% than those of the explicit calculation. 

V. CONCLUSIONS 

A large scale investigation of the three-body dispersion coefficients for a set of 14 atoms has 
been completed. The results presented for the set of small systems, namely H, He and Li should 
be regarded as the existing calculational standard. 

The current set of values should also be adopted as the reference standard to be used for heavier 
systems. Even though more approximate methods of solving the Schrodinger equation are adopted, 
the semi-empirical Hamiltonian adopted for the alkali atoms and the alkaline earth atoms reduces 
the computational complexity to such an extent that uncertainties in the solution of the resulting 
Schrodinger equation are largely eliminated. When comparisons of polarizabilities and two-body 
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TABLE VII: The three-body dispersion coefficients (in atomic units), for systems containing two neon 
atoms and two argon atoms. The sources for the /-values distributions are the same as those in Table IVl 
Results in the Midzuno-Kihara (MK) approximation 



35 



use the present oscillator strength distributions 



to compute the underlying ay and Cq needed as input. The numbers in the square brackets denote powers 
of ten. 



Systems 


Present 


MK 


^112 


•^211 


•^122 


•^221 


^113 


•^311 










X-Ne-Ne 










H 


12.64 


13.01 


35.61 


49.82 


100.3 


140.4 


200.5 


467.3 


Li 


111.9 


113.6 


315.2 


1907 


887.9 


5370 


1768 


6.502[4] 

L J 


Na 


126.2 


130.2 


355.6 


2567 


1002 


6638 


1994 


8.629[4] 


K 


191.6 


200.9 


539.9 


4803 


1521 


1.353[4] 

L J 


3028 


2.268[5] 

L J 


Rb 


215.1 


228.0 


606.0 


5785 


1707 


1.630[4] 


3399 


2.893[5] 


Be 


62.79 


64.37 


176.9 


674.8 


498.3 


1901 


994.2 


1.060[4] 


Mg 


102.5 


105.6 


288.8 


1491 


813.5 


4199 


1622 


3.156[4] 

L J 


Ca 


175.6 


183.1 


494.6 


175.6 


1393 


1.073 [4] 


2776 


1.157[51 

L J 


Sr 


212.2 


223.0 


598.0 


5111 
X-Ar-Ar 


1685 


1.440 [4] 


3356 


1.754[5] 


H 


174.0 


177.9 


846.6 


658.1 


4127 


3217 


8470 


6056 


Li 


1818 


1840 


8704 


2.955 [4] 


4.170[4] 


1.420[5] 


8.804[4] 


9.836[5] 


Na 


2022 


2068 


9690 


3.675 [4] 


4.647[4] 


1.766 [5] 


9.793[4] 


1.314[6] 


K 


3041 


3169 


1.459 [4] 


7.636[4] 


7.006[4] 


3.662[5] 


1.474[5] 


3.532[6] 


Rb 


3378 


3557 


1.622 [4] 


9.220[4] 


7.798[4] 


4.208[5] 


1.638[5] 


4.525 [6] 


Be 


945.6 


962.3 


4557 


9714 


2.198[4] 


4.704[4] 


4.585 [4] 


1.476[5] 


Mg 


1568 


1600 


7544 


2. 220 [4] 


3.633[4] 


1.071 [5] 


7.599 [4] 


4.548 [5] 


Ca 


2719 


2812 


1.307[4] 


5.890[4] 


6. 290 [4] 


2.832[5] 


1.318[5] 


1.732 [6] 


Sr 


3279 


3415 


1.577[4] 


7.970 [4] 


7.590[4] 


3.829[5] 


1.590[5] 


2.655[6] 



dispersion coefficients have been made the CICP approach has 



34 



typically been within 1% of the 



54]- 



most sophisticated approaches to atomic structure available 

Besides an improvement in accuracy, the scope of the present work exceeds that of any previous 
work by at least an order of magnitude. For example, Huang and Sun [45J, and Patil and Tang 48]. 



presented extensive results for the alkali atoms. However, the present calculations also encompass 
the rare gases and the alkaline-earth atoms. 
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TABLE VIII: The three-body dispersion coefficients (in atomic units), for systems containing two krypton or 
two xenon atoms. The sources for the /-values distributions are the same as those in Table [Vj Coefficients 
in the Midzuno-Kihara (MK) approximation 35, use the present oscillator strength distributions to 
compute the underlying ad and Cg needed as input. All values are in atomic units. The numbers in the 
square brackets denote powers of ten. 



Systems 


Zm 


^ 112 


■^21 1 


' J 122 


"221 


^ 113 


^311 




Present 


MK 






















X-Kr-Kr 










H 


370.1 


380 


2485 


1385 


1 352 [41 


8409 


2 760 [41 


1 268 [41 


Li 


4081 


4136 


2.407[4] 


6 612 [41 


1 420 [51 


3 870f5l 


3 044 [51 


2 151[6l 


Na 


4521 


4628 


2.671 [4] 


8.216[4] 


1.579[5] 


4.819[5] 


3.372[5] 


2.881[6] 


K 


6795 


7.086 


4.017[4] 


1.700 [5] 


2.377[5] 


1.006[6] 


5.070[5] 


7.800[6] 


Rb 


7526 


7934 


4.454[4] 


2.055[5] 


2.639[5] 


1.215[6] 


5.616[5] 


1.000[7] 


Be 


2066 


2109 


1.232[4] 


2.092 [4] 


7.354[4] 


1.256[5] 


1.539 [5] 


3.147[5] 




3445 


3523 


2 049 [41 


4.831 [4] 


1 220 [51 


2 885 [51 


2 567f5l 


9 796 [51 


Ca 


6015 


6229 


3 569 [41 


1 299 [51 


2.120[5] 


7.714[5l 


4.484[5] 


3.774[6] 


Sr 


7253 


7563 


4.304[4] 


1.763 [5] 
X-Xe-Xe 


2.557[5] 


1.046 [6] 


5.407[5] 


5.807[6] 


H 


871.4 


899.3 


7139 


3217 


5.862 [4] 


2.650[4] 


1.139[5] 


2.930[4] 


Li 


1.034[4] 


1.049 [4] 


8.232 [4] 


1.610[5] 


6.558[5] 


1.294[6] 


1.326[6] 


5.255[6] 


Na 


1.140 [4] 


1.168[4] 


9.092 [4] 


2.013[5] 


7.257[5] 


1.616[6] 


1.464[6] 


7.059 [6] 


K 


1.713[4] 


1.787 [4] 


1.366[5] 


4.528[5] 


1.091 [6] 


3.405 [6] 


2. 200 [6] 


4.258 [7] 


Rb 


1.891 [4] 


1.994 [4] 


1.510[5] 


5.157[5] 


1.208 [6] 


4.120[6] 


2.430[6] 


2.483[7] 


Be 


5042 


5162 


4.072[4] 


5.009[4] 


3.293[5] 


4.076 [5] 


6.531[5] 


7.439[5] 


Mg 


8476 


8690 


6.824[4] 


1.172[5] 


5.500[5] 


9.489 [5] 


1.095[6] 


2.345[6] 


Ca 


1.494 [4] 


1.549 [4] 


1.199[4] 


3.210[5] 


9.625[5] 


2.580[6] 


1.962[6] 


9.177[6] 


Sr 


1.802 [4] 


1.881 [4] 


1.445 [5] 


4.377[5] 


1.161 [6] 


3.512[6] 


2.323[6] 


1.419[7] 
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